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Abstract 

We investigate the short-, medium-, and long-term time dependence of wave packets in the 
infinite square well. In addition to emphasizing the appearance of wave packet revivals, i.e., 
situations where a spreading wave packet reforms with close to its initial shape and width, we 
also examine in detail the approach to the collapsed phase where the position-space probability 
density is almost uniformly spread over the well. We focus on visualizing these phenomena in both 
position- and momentum-space as well as by following the time-dependent expectation values of 
and uncertainties in position and momentum. We discuss the time scales for wave packet collapse, 
using both an autocorrelation function analysis, as well as focusing on expectation values and find 
two relevant time scales which describe different aspects of the decay phase. In an Appendix, 
we briefly discuss wave packet revival and collapse in a more general, one-dimensional power-law 
potential given by Vi^k-^{x) = VQ\x/a\^ which interpolates between the case of the harmonic oscillator 
{k = 2) and the infinite well {k = oo). 



I. Introduction 



The study of wave packet solutions of the Schrodinger equation in model potentials 
can illuminate many aspects of wave mechanics, both semi-classical features which have 
an obvious classical analog, as well as purely quantum mechanical effects. Since the first 
numerical studies of wave packet motion [l| appeared in this journal, there have been a 
large number of papers illustrating such phenomena as wave packet interactions with square 



barriers or wells j^, multi-well systems iJl], multi-dimensional scattering j^, and in systems 
of relevance to solid-state physics [5j. Popular simulation packages lc| now allow students 
to easily visualize the time-evolution of quantum states in a number of systems (as opposed 
to the more traditional static images of single stationary state solutions) and allow them to 
systematically vary the values of important physical parameters to study the dependence 
on such quantities as the particle mass, incident energy, and barrier/well properties such as 
height/depth and width. 

Such studies have become increasingly relevant as pedagogical tools as wave packet prop- 
agation in more realistic quantum mechanical bound state systems has been probed exper- 
imentally, especially the behavior of Coulomb wavepackets on circular |7| or elliptical ^ 
orbits which are accessible in Rydberg atom states jsf. One of the most intriguing aspects 
of such systems is that initially localized wave packets, which have a short-term time evo- 
lution which exhibits simple classical behavior, will spread significantly after several orbits, 
entering a so-called collapsed phase, only to reform later in the form of a quantum revival 
in which the spreading reverses itself and the wave packet relocalizes. 

One of the early papers to point out the possibility of such collapse/revival behavior in 
the context of the infinite well was by Segre and Sullivan They studied the motion of 
bound state wave packets in such a system using a simple "sum over energy eigenstates" 
method. These authors noted, in passing only, many of the same aspects observed in the 
Coulomb system, namely the spatial spreading of the wave packet, continuing until the 
probability density is almost uniformly spread over the entire well, and then reforming into 
something like the t = configuration. 

More recently, Bluhm, Kostelecky, and Porter have examined the time evolution 
and revival structure of wave packets in more generic one-dimensional systems (examining 
the harmonic oscillator and infinite well in particular as special cases) while Aronstein and 



Stroud |3l have focused on fractional revivals [3| in the infinite well, also presenting visu- 
alizations of position-space probability densities. These studies have made extensive use of 
the autocorrelation function introduced by Nauenberg [3], defined by 



C{t) = I ip*{x,t)'^{x,0)dx 

<f)*{p,t)<j){p,0)dp (1) 
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where the final form is useful if the bound state wave packet is written in terms of energy 
eigenstates, Un{x), via 

oo 

=^a„M„(x)e-*^"*/'^ (2) 

n=0 

The quantity C{t), which measures the overlap of the initial state wave packet (in either 
position- or momentum-space) with the state at later times, can be used in its last form to 
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times 
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"ormation about the classical period of motion, as well as revival and superrevival 



using information on the n-dependence of the En- 
While the autocorrelation function is a powerful tool for such formal and general analyses 
iJ ]. plots of the time-dependence of C{t) are rather distant from the intuitive picture of 
initially semi-classical particle motion in the well, accompanied by wave packet spreading, 
leading to a collapsed state. Plots of position-space probability densities, such as those in 
Ref. [IJ] for isolated times, contour plots {\ip{x,t)\'^ versu^(x, t)) for all t up to a revival 
period (as in the 'quantum carpet' visualization in Ref. IIgJ), or even in the form of more 
dynamic animations, all can give valuable information on the behavior of the quantum state 
in position-space as a function of time. 

A much more direct connection, however, to an intuitive classical description involving 
particle trajectories can be obtained by evaluating the time-dependent expectation values of 
standard variables such as position and momentum, {x)t, {p)t, which can then be compared 
to familiar classical results, illustrating the expected quasi-classical behavior over the first 
few periods, but also demonstrating the results of purely quantum effects during the collapse 
and various revival phases. Since the collapse is due to the spreading of the wave packet, it 
is also natural to 'track' the behavior of the spreads or uncertainties in the same variables, 
Axt, Apt, as functions of time, examining as well how these behave during the semi-classical, 



collapse, and revival phases. Taken together, these can offer a natural visualization and 
complementary description of the quantum state over all of the natural time scales in the 
problem, including the semi-classical and collapsed phases which have not been examined in 
as much detail as the important revival structures. While the results we present here will be 
evaluated numerically, our approach is the same as that taken by Styer [l7| in understanding 
wave packet propagation through expectation values and uncertainties. 

In this note, we begin by very briefly reviewing the autocorrelation approach to the re- 
vival structure of wave packet propagation in the infinite well in Section IIA. We illustrate 
this with calculations, not only of the standard C{t) over the revival period, but of a related 
correlation function which provides information on the reformation of the wave packet at 
half-revival periods as well as providing additional illustrative evidence for fractional re- 
vivals. Probability density plots, in both position- and momentum-space, are also provided 
to illustrate the connections to the autocorrelation function approach as we review and ex- 
tend the results in Refs. Q and Q, but also for later comparison to results obtained from 
expectation values. In Section IIB we focus on the short-term evolution of wave packets, 
visualizing through expectation values and uncertainties in both x and p, the quasi-classical 
propagation and the quantum spreading leading to collapse. In Section IIC, we repeat this 
analysis, but instead focus on the entire revival period, providing more graphic evidence of 
the reformation at the half-revival time. In contrast to earlier works, in Section IID, we focus 
on the medium-term time dependence, emphasizing the approach to the collapsed phase. We 
study the decay in the degree of correlation between the initial packet and later states (as 
measured by |C(t)|) as well as the approach to the quasi- uniform or 'fiat' probability density 
phase discussed by Segre and Sullivan jl^. We find that the time scale for the initial loss 
of coherence as measured by \C(t)\ is simply related to the free particle spreading time, to, 
while that which describes the approach to the 'fiat' or truly collapsed state, tg^^^ (as given 
by the approach of {x)t, {p)t, and most obviously Axt to their 'fiat' values), are, in fact, 
different and we examine both in detail. In this regard, we find that an approach involving 
expectation values and uncertainties provides useful information which is not encoded so 
obviously in the autocorrelation approach. 

Finally, in an Appendix, we consider wave packet revivals in a more general, one- 
dimensional power-law potential given by V(^k){x) = VQ\x/a\^ as this form interpolates 
smoothly between the case of the harmonic oscillator {k = 2) and the infinite well {k = oo). 
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Usin g si mple WKB methods for the estimation of the energy eigenvalues (as was done in 
Ref. jisl) and the methods in Ref. we can explicitly exhibit the scaling of the revival 
time with k which gives the familiar result for the infinite well and which also exhibits the 
appropriate divergence [Trev — > oo as A; — » 2) as one approaches the harmonic oscillator 
case. For this same general system, we also evaluate the revival times for half-wells where 
one introduces an infinite wall at the origin and restricts the particle in V(fc)(x) to have 
X > 0. Gea-Banacloche has recently examined the case of a 'quantum bouncing ball' in 
such a potential {V{z) = mgz for 2; > so that k = 1) and we find that our general result 
reproduces that special case. 

II. Visualizing wavepacket revivals and collapse 



While the revival structure of wave packets in the infinite well has been discussed in 



detail in Refs. 



ill and 



12| , we will briefly recall the most important points of an autocorre- 



lation function analysis of this system, extending these results slightly and providing some 
complementary visualizations of related topics. 



A. Autocorrelation function analysis 



For definiteness, we consider the one-dimensional infinite square well system, with walls 
at a; = 0, L, for which the most general time-dependent solution of the Schrodinger equation 
(since this system does not allow for unbound or continuum states) is given by the infinite 
discrete summation 

00 

^(x,t) = ^a„M„(a;)e-^^"*/" (3) 

n=l 

where the (normalized) eigenstates and energy eigenvalues are, of course, 

/A /2" . fmTx\ , ^ 2 ^^^^ /A\ 

Unix) = y ;^ sm (^-^ j and = n ^— ^ (4) 

Recall that the normalization of the position- and momentum-space wavefunctions are re- 
lated via 

l = $^|an|'= / mX,t)\^dx= / \(P{p,t)\''dp=l (5) 
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Because we are interested in evaluating expectation values of many quantities for somewhat 
detailed comparison, we will be careful to ensure that any wave packet is appropriately 
normalized. 

Because the quantized energy eigenvalues are integral multiples of a common value, it is 
easy to see that any such solution will exhibit revivals or reformations, where ip{x,t + T) = 
ip{x, t), when 

e-'^"^/^ = +l or r^^^rnL^ .g. 

n n 

which will be called the revival time. As shown in Refs. and the wave packet will 
also reform itself at half this time as well, but in a possibly different location. We see this 
by noting that 



ip{L - x,t + T /2) = anUnjL - x)e~ 



{t+T/2)/h 



oo 



= anUn{x)e-''^-'/^ [- cos(n7r) cosin^T:)] (7) 

n=l 

This connection can be written in the form 

^{L-x,t + T/2) = -^{x,t) =^ ^{x,t + T/2) = -ilj{L-x,t) (8) 

or 

\ij{x,t + T/2)\^ = \ij{L-x,t)f (9) 

so that at half the revival time later, any initial wavepacket will reform itself (same shape, 
width, etc.), but at a location mirrored about the center of the well. 
The momentum behavior of the system will be given by 

1 

0(p, t) = -^== J V^(x, t) e^f^/^ dx (10) 

at all times and we can use this connection to understand the behavior of the wave packet 
in momentum space at T/2. We note that 

1 '•^ 



0(p,t + T/2) = -j== I i){x,t + T/2)e''P'''^dx 
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ij{L-x,t) e'^'^'^dx 

^{y,t) e-'^y'^dy 



L 







(11) 



so that 

|0(p,t + T/2)p = |0(-p,t)r (12) 

and half a revival time later the initial momentum profile is also reproduced, except flipped 
in sign {p —>■ —p), so that the particle is moving in the 'other direction'. These statements 
are true of any general time-dependent solution of the form in Eqn. (jS)) for the infinite square 
well. 

For a typical wave packet solution, the expansion coefficients are such that they are 
sharply peaked about some large value of n, say uq » 1, so that the packet can be 
characterized by a momentum value po ~ n^Tih/ L. The period of the classical motion will 
be given by r = 2L/fo and if we associate vq = pq/itl we find that 



_ 2L _ 2L'^m _ 1 
Po/m riQTch 2nQ 



AmL 
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Pn/iii un'iuL zilq iiii 2rin 
(See also Ref. [ill for a more general discussion of the classical period.) We thus expect the 
wave packet to undergo many classical periods before reforming itself in a revival and note 
that there are at least two important time scales in the problem, namely T >> r. 

As an explicit example of a bound-state wave packet, we will use a quasi-Gaussian form, 
with expansion coefficients given by 



V L 
where 

rnih , nQTih 

Pn = -j^ = (1^) 

define the central momentum, corresponding to an initial speed to the right oi vq ^ po/m. 
The initial position of the wave packet is then given by xq. (This is similar to the form used 
in Ref. Q-) 

This form gives an initial wave packet which is almost Gaussian and hence very similar 
in most regards to the standard analytic result seen in many textbooks where explicit free- 
particle wave packets are constructed. (See, for example, Refs. [20 1, 21 1.) For example, 
with this normalization, the expansion coefficients are already almost normalized (but we 
explicitly ensure the proper normalization in all of our numerical calculations) and the initial 
uncertainties or spreads in the wave packet are given (to an excellent approximation) by 



provided Axq << L. We recall that the corresponding free particle wave packet with 
Gaussian shape given by the continuous momentum weighting 



= J_^^-aHp~Por/2^-ipxo/h ^^7^ 



spreads with time with a spatial uncertainty given analytically by 



Axt = Axo^/l + ( where to = mha'^ = ^"^(^^o) (ig) 

This spreading time, to, defines a third time scale in the problem which we can write in the 
following form 

^.„.(^)^ (19) 

for future comparison. For our explicit calculations, we choose the numerical values 2m = 
h = L = 1 to define the model system and for the specific wave packet solution we use 

Axq = 0.05 — > a = -^-^ — y Ap = 10 (20) 

10^2 ^ ^ 

with 

xo = 0.5 and no = 400 (21) 

With these values we have T/t = 2nQ = 800 and to/r = tt so that the wave packet will 
spread significantly on the time scale of a few classical periods. 

Using this explicit wave packet, we can evaluate the autocorrelation function given by 
Eqn. and we plot |C(t)| versus t over half a revival period in Fig. 1. (The plot is 
symmetric about the half-revival time, T/2.) The initial decrease in correlation for t > 
is seen as the 'decay' in |C(t)|, while the structures in C{t) due to many of the fractional 
revivals jl^ are also apparent. In order to show evidence for the reformation at T/2 
mentioned above, we also show in Fig. 1 a plot of a related correlation function defined by 

/+00 
^*{L-x,t)tlj{x,0)dx 
■oo 

hoo 



^+ipL/h ^22) 



which measures the correlation between the initial wave packet and later states which are 
'fiipped' about the center of the well or, equivalently, reappear with opposite momentum 
profiles. Using this form, the reformation at T/2 is clearly seen (|C(T/2)| = 1) as well as 
additional evidence for fractional revivals. 
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Also using this specific wave packet as an example, we illustrate the position-space prob- 
ability densities at several times in Fig. 2. The initial, centered wave packet at t = (solid 
curve) can be seen later at several obvious fractional revival times. The packet at t = r/8 
(dashed curve) can also be seen one-half revival time later, reformed on the other side of the 
well as dictated by Eqn. Q. The position-space probability density at a more 'random' later 
time (if: = 124r) illustrates the observation of Ref. where the position-space probability 
has become almost equally distributed over the entire well; the arrow indicates the magni- 
tude of a truly uniform probability distribution given by Pcl{x) = 1/L (corresponding to 
the numerical values used here.) 

In order to extract information on the momentum-space wavef unctions, we can Fourier 
transform Eqn. Q, but we can also formally invert ilj{x,t) to write 

oo 

0(p,t) = ^a„0„(p)e-^^"*/^ (23) 

n=l 

where the momentum-space eigenstates corresponding to the Un{x) are given by 




so that 

/ + 00 
p\<p{p,t)\^dp (25) 
'OO 

and similarly for other powers of momentum, which will be useful below. 

We plot, in Fig. 3, the momentum-space probability densities corresponding to the times 
in Fig. 2 and note that the wave packet is characterized by a reversal in momentum at T/2 
as in Eqn. (fT^ . The initial value of the momentum spread {Ap = 10) is also recovered at 
the half-revival, but the uncertainty at many later times is governed not by the intrinsic 
width of the initial (f){p, 0) peak, but rather by the 'distance' between the two sharp and 
well-isolated momentum peaks at approximately ±po which gives Ap ^ +pq for much of the 
collapsed phase. 

While these examples are somewhat illustrative of the time-dependence of the expectation 
values and uncertainties in position and momentum variables, we now turn our attention to 
more detailed examinations of the behavior of {x)t, Axt, {p)t, and Ap^ over three important 
time scales. 



B. Expectation value analysis: short-term, classical behavior 

9 



In many situations, the hierarchy of times given by r < to << ^ will hold and we will 
examine the time-dependence of the wave packet propagation over the first few classical 
periods. Using our standard parameter set, we calculate the expectation values and spread 
in position over a time scale equal to the first ten classical periods, (0, lOr) and plot the 
results in Fig. 4. The periodic and almost classical nature of the time-dependence of {x)t is 
clear (compare to Fig. 5 which gives the purely classical result for x{t),v{t) over the same 
period), but obvious differences are also present. The amplitude of the motion, operationally 
defined by half the total "left-to-right" distance in {x)t (which for a purely classical particle 
would be constant and equal to L/2) clearly decreases with time in the quantum case and the 
reason for this is clear from the corresponding plot of the time- dependent spatial uncertainty, 
Axf. The wave packet disperses in time (away from its initial Axq = 0.05 value, given by 
the horizontal dashed line) and its width increases, initially at least, in a manner completely 
consistent with a free particle Gaussian wave packet, as the dotted curve which forms the 
upper envelope for the solid data is given by Eqn. (fT8|) . (In this case, the importance 
of including the exact quantum time- dependence and not ignoring the dispersion is 
obvious.) The sharp 'dips' in Axt correspond to 'collisions' of the wave packet with the 
walls and have been recently studied in the context of a single 'bounce' 23|. As the wave 
packet gets wider, the value of {x)t at the classical collision times decreases (at the right 
wall) or increases (at the left wall) as the spreading wave packet finds it more difficult to 
approach the walls. Despite this obvious quantum effect, the short-term time dependence 
of the system is relatively classical, just as in the Coulomb case. 

Turning attention now to the momentum-space behavior, we can make use of Eqn. ()2H|1 
to evaluate {p)t and Apt and we plot these values in Fig. 6 over the first ten classical periods. 
The correspondence between the sharp classical 'jumps' in the v{t) versus t plot in Fig. 5 and 
the t-dependent quantum expectation value of p in Fig. 6 is clear, but these changes become 
increasingly 'softer' as the wave packet broadens. The plot of Am versus t is consistent with 
expectations from the 'single-bounce' case considered in Ref. j2^. During periods when 
the packet is not involved in collisions, the momentum spread is initially simply given by 
the intrinsic width of (f){p), ApQ, but increases dramatically during a collision with a wall. 
At such points, where the momentum distribution is shifting from a single large peak at 
positive p = +Pq to a similar one at p = —po, the spread in momentum values is dominated 
by the 'distance' between the two peaks. For example, a simple model for the momentum 
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distribution at the collision time might be 

Pcl{p) = \ [6{p - Po) + 6{p + po)] (26) 

which gives {p) = and Ap = +po. This behavior is clearly seen during collisions for times 
satisfying t < to. The increasing spread in the position-space wave packet (approaching 
Ap ~ +Po) clearly implies that the momentum distribution comes to include significant 
components with both p ^ ±po for later times as well, giving the same large value of Apt 
over much of the collapsed phase (as seen earlier in the numerical values in Fig. 3.) 

C. Expectation value analysis: long-term revival behavior 

Turning now to the long-time dependence of the wave packet, we plot, in Fig. 7, {x)t 
versus t over one entire revival time (T = 800r.) In the top plot of Fig. 7, we show {x)t 
for the wave packet at times given by t = {n + l/8)r and {n + 5/8)r as the solid and 
dashed curves respectively. (The t = nr data yield the horizontal line at {x)t = 0.5L.) The 
(n+ l/8)r data (solid curve) corresponds to those times at which the classical particle would 
be one-half of its way from the center to the right-hand side and classically would always be 
located at {x)t = 0.75L, as indicated by the horizontal dashed line: the (n+5/8)r case would 
similarly give {x)t = 0.25L in the classical limit. While the wave packet data are initially 
consistent with these classical results, they quickly flatten during the collapsed phase so that 
the expectation value of position at even these times is consistent with (x) = 0.5L. We note 
that at t ~ 800r = T the expectation values again match the classical results, which is our 
first indication of the expected revival or reformation from this expectation value approach. 

Similarly, we note that at t ~ 400r = T/2, we find that the expectation values for position 
have reversed themselves, consistent with Eqn. (jH)) for the half-revival. In the lower half of 
Fig. 7, we show results for {x)t for all times in the range t = (116r, 164r) to illustrate the 
fine structure apparent in the time-dependence and how the wave packet oscillates, with 
very small excursions, around {x)t = 0.5L during most of the collapsed phase. 

In order to further probe the manner in which the probability distribution 'fiattens' during 
the collapse and reforms during a revival, we plot in Fig. 8 the uncertainty Axt versus time 
over the same time periods as in Fig. 7: we show values for times given hj t = nr only. To 
compare the results to a putative fiat probability distribution, we note that this case would 
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be given by 
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(27) 
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and would yield expectation values 
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We note that Axj grows (consistent with Eqn. (fTH|l ). but rapidly flattens out during much of 
the collapsed phase at a value which is indeed consistent with Axg^^^ = 0.288L over much of 
the entire revival period. (At the bottom of Fig 8, we show the behavior of Aa;j for all times 
in the interval (116r, 164r) indicating the oscillations around the expected 'flat' result.) 

The revivals at T = 800r and T/2 = 400r are also apparent here as the uncertainty 
returns to its initial value, indicating a reformation of the initial wave packet, perhaps at a 
different location. We also note, however, the additional feature of quasi-revivals at t = 200r 
and 600r (T/4,3T/4) where the wave packet returns to its initial width. To illustrate the 
nature of the wave packet at these times, we have already shown in Fig. 2 the form of 
the packet near t ^ 200r. The wave packet at t = 200r does not precisely return to its 
initial state (the real part of ip{x,T/4) is larger than ip{x,0), while the imaginary part of 
iIj{x,T/4:) vanishes), but does so in a way which reproduces Axt=T/4 = Axq. We can see 
that this is a very special CclSG, clS the t = t/8 wave packet is obviously split into two pieces 
at t = r/8 + T/4 as it is 'half way' to being reformed at the half-revival time. We can 
also note that at t = lOOr, 300r, 500r, and 700r that there are 'anti-revivals' in that the 
wave packets are spread somewhat more than during the rest of 'flat phase' of the collapsed 
phase. The position-space wavefunction at t = lOOr was also shown in Fig. 2 to illustrate 
this 'anti-correlation'. 

We can also visualize the long-time behavior of the momentum-space variables by plotting 
{p)t and Apt versus t over one revival time in Fig. 9 (again, for values of t = nr.) We note 
that the momentum spread quickly saturates at the 'flat' value of Ap = +Pq, except when 
it returns for revivals at t = T/2 and T. The expectation value of momentum reverses sign 
at t = T/2, as expected from Eqn. (fT^ . while returning to its original value only at the 
standard revival when t = T. 

While we have explicitly considered only the quasi-Gaussian wavepacket with specific 
numerical values here, many of the same effects are seen quite generally. Changes in most 
of the parameters simply change the appropriate time scales involved. Wave packets with 
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different initial positions, xq 7^ L/2, do not exhibit the special quasi-revivals at t = T/4, 3T/4 
because they lack the required additional symmetry. Different functional forms for the 
expansion coefficients (the discrete versions of Lorentzian momentum weightings or initially 
flat spatial wave packets, for example) all exhibit similar behavior to that discussed here. 

D. Expectation value analysis: medium-term, collapse behavior 

One aspect of the behavior of wave packets in the infinite well which has not been ex- 
amined in much detail as has the pattern of revivals (fractional or otherwise) is the decay 
phase, the time over which the spreading of the wave packet leads to an increasing lack of 
coherence (as measured by C(t)) and because of which the position-space probability den- 
sity approaches the 'fiat' distribution observed by Segre and Sullivan One measure of 
a decay time scale associated with this spreading is given by the scaling behavior of |C(t)| 
at integral values of the classical period, namely t = nr. For the oscillator, for example, 
we have \C{t = nr)| = 1 since the wave packet motion reforms exactly with the classical 
oscillation frequency. For the infinite well, we find (after a large number of numerical trials) 
that the autocorrelation function for such integral times is initially given by 

|(:7(t = nr) I ~e-(""/^^)' for nr < Tc (29) 

where the collapse time, Tq, is given by 

and where An is the dispersion in the n-distribution for the Gaussian wave packet. This 
result may not be unexpected since the decrease in coherence measured by C{t) is directly 
tied to the spreading of the wavepacket and to is therefore the natural time scale which arises 
in an autocorrelation analysis of this type. 

Examination of Fig. 8, however, suggests that there is a second scale involved in the decay 
process, namely that associated with the time it takes for the wave packet to approach the 
spatially fiat probability distribution. In order to better visualize the approach to the 'fiat' 
distribution, we plot Axt versus t in the range (0,T/8 = lOOr) for three different values of 
Axo in Fig. 10. We note that in each case the uncertainties initially follow the free-particle 
Gaussian result in Eqn. (fTHj). but then level off to the same 'fiat' value of Axgg^^ = L/\/T2. 
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We can estimate the time it takes to reach the constant value during the collapse phase by 
equating 

which gives 

which defines still another time scale, typically an intermediate scale between to and T. 
Numerical studies, such as in Fig. 10, indicate that the position-space uncertainty saturates 
at its 'flat' value at a time given roughly by Sigg^^.. The three cases shown in Fig. 10 then 
illustrate that this flattening time does indeed scale with Axq, in contrast to the coherence 
decay time scale Tc which is proportional to Axq. We can see the same flattening time scale 
(and its proportionality to Axq) by plotting the medium-scale time dependence of {x)t, {p)t, 
and Apt in Fig. 11, all of which approach the appropriate 'flat' value with a time scale of 
roughly St^g^-^.. It is not immediately apparent how information on this important time scale 
is readily obtainable from an autocorrelation analysis (for example, as shown in Fig. 1), 
providing more evidence for the usefulness of this more intuitive approach. 

E. Conclusions 

We have studied the time-dependence of wave packets in the inflnite well, their collapse 
and revival structure, and the time scales appropriate for each phase by examination of the 
expectation values and uncertainties in the familiar variables of one-dimensional quantum 
mechanics, position and momentum. The visualizations provided here complement more 
formal derivations of revival structures studied using autocorrelation function analyses. In 
the case of the approach to the collapsed phase of the wave packet approaching a 'flattened' 
state, this type of expectation value analysis very naturally yields information on the scaling 
properties of the natural time scale for the collapse, in contrast to an autocorrelation function 
approach where such information is not so immediately and intuitively obvious. 
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Appendix 



Motivated by the simple revival structure observed for the harmonic oscillator and infinite 
wen pote... . Ref. A we p.e.n. .e.e a bnef a„.... of .Ke ™, ..e to. .Ke .ene.al. 
one-dimensional power-law potential 

• k 



V^k){x) = Vo 



(33) 



which interpolates between these two special cases for /c = 2 (oscillator) and k = oo (infinite 
square well). While there is really only one independent dimensionful parameter in this 
potential, namely Vo/a^ = a, we write it in this form to emphasize the k (yo limit where 
we obtain an infinite well (of width 2a, compared to our earlier analysis.) A simple variation 
is the 'half power-law potential given by 



I V(fc)(x) forx>0 
10 for a; < 



(34) 



which is then appropriate for the analysis of the 'quantum bouncing ball' discussed recently 
in Ref. when k = 1 and Vo/a 



mg. 



In order to make use of an autocorrelation analysis as in Ref. we require the quantized 

(k) 

energy eigenvalues, En ■ Since we are dealing with situations in which n » 1, it will be 
a good approximation to utilize the WKB energy quantization condition to evaluate the 
Eit^ . The standard version (which assumes wave-function matching at linear walls via Airy 
functions) is given by 



+X0 



XQ 



2m 



Ei'^ - V( 



(k) 



with n 



0, 1, 2, ... and where 



X) 



Er, 



dx 



l/k 



n H — \ hn 
2 



(35) 



(36) 



are the classical turning points. This expression is known to be exact for the harmonic 
oscillator. When the boundary conditions at a given turning point are imposed by infinite 
walls where the wavefunction must actually vanish, the appropriate matching coefficients 
are Cl,Cr = 1/2 instead of the Cl,Cr = 1/4 for linear walls 2J|, so that for the infinite 
square well the right hand side is replaced by (n + 1/2) —>■ {n + 1) which also gives the exact 
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answer for this case. (The combination of the two to describe the 'half oscillator' where 
V{x) = for x > and an infinite wall at the origin is then also given exactly 

with these two cases.) Using the standard WKB form for all /c < oo, we find the energy 
quantization condition 



(n+ 1/2) 



r(i/fc + 3/2) 

2av^'" r(l/A; + l)r(3/2) 



1/k 



2k/(k+2) 



(37) 



(where {n + 1/2) {n + 1) for k oo.) This reproduces the exact oscillator and infinite 
well examples in the k = 2 and k ^ oo limits. The classical periods for a given value of k 
and n will then be given by 



-ik) 



2Tih, , , f2 + k 
(n + 1/2) 



(k) 



2k 



T{k, n) 



(38) 



which reduces to the familiar oscillator and infinite well results. Using the result of Ref. 
for the revival time, namely 



ink 



(PEn/dn'^ El 
we find that the revival times are given by the simple expression 

A; + 2 



T 



T{k,n) 



k-2 



2nr(A;, n) 



(39) 



(40) 



This can be confirmed by explicit evaluation of |C(t)|, \C{t)\ over one revival time for any 
value of A; > where one finds plots which are very similar to Fig. 1, except that the revivals 
are no longer exact as they are in the case of the infinite well. This form explicitly exhibits 
the divergence we expect in the oscillator case (when k 2) where the wave packets are 

nn 

exactly periodic [25||, |26l |. 

We have also performed a more numerical analysis using the quasi-Gaussian wave packet 
to examine the short-time decay in correlations as measured by the decrease in the magnitude 
of |C(t)|. When we evaluate the autocorrelation function at integral values of the classical 
period, we again find that 



where 



c 



\C{t 
k + 2 



nT)\ ~ e 



k-2 



-{nr/Tcf 



1 



T{k,n) 



27rAn2 27rAn2 



(41) 



(42) 



2nT(k, n) 

which reproduces the result for the infinite well we have considered more explicitly, as well 
as exhibiting the appropriate divergence as A; ^ 2 as we approach the exactly periodic 
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oscillator case. This result also holds for other, non-Gaussian wave packet expansions, 
with the numerical factor 2tt replaced by different numerical coefficients of the same order. 
Two open questions which would be interesting to explore further in this context of the 
general well are whether during the collapsed phase the position-space probability naturally 
approaches a semi-classical probability density (as in Ref. over much of the collapsed 
phase and, if so, how the 'flattening' time to approach this semi-classical distribution scales 
with k and other parameters, i.e., what is the generalization of Eqn. ()32|) . 

The case of the 'half general-power law potential given by Eqn. ()34|) is easily analyzed 
in the same manner. The WKB energy quantization condition is applied at the classical 
turning points x = +Xo and a; = with the result 



(n + 3/4)^=^0^ 



hn ^^yk r(l/fc + 3/2) n2/c/(A.+2) 



(43) 



aV2^ ' r(l/A; + l)r(3/2)_ 

with the factor (n + 3/4) {n+ 1) for k —>■ oo. (Recall the discussion after Eqn. (jHBj) about 
the appropriate matching coefficients for infinite wall boundaries.) The classical period is 
then given by Eqn. ^ with E^n^ E^n^ and {n + 1/2) ^ (n + 3/4). We find that the 
same relationship between revival time and classical period as in Eqn. (jin|l holds and we 
note that our general result then reproduces the revival time found in the special case of the 
'quantum bouncer' considered in Ref. 19| when k = 1 and V^/a = mg. 



Gea-Banacloche 



19| has also discussed a collapse time for the case of the 'quantum 



bouncing ball' (the k = 1 case of the 'half-well') and we note in passing that his expression 
is close in spirit to the 'flattening' time discussed for the infinite well as his expression for 
the collapse time scales as 1/An as does tf^^^, and not, for example, like Tc oc 1/An^ for 
the decay time obtained by the scaling properties of the autocorrelation function. 
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Figure Captions 



Fig. 1. Plot of the autocorrelation function, |C(t)|, (in Eqn. (J])) for the infinite square well. 
The modulus of C{t) is evaluated at integral values of t/r over the first half revival 
time; the curves are symmetric about t = T/2 = 400r. The vertical dashed lines 
indicate values of some of the possible fractional revivals, namely {p/q)T for integral 
values of p < q. The standard wave packet parameters in Eqns. (HH), (plj) . and fETj) 
are used. The related correlation function, C{t), defined in Eqn. ((221), is also plotted 
and shows much of the same structure, including the 'anti-revival' at t = T/2 = 400r 
where the wave packet reforms, but with opposite momentum values. 

Fig. 2. Position-space probability distributions, |?/'(a;,t)p versus x, for the quasi- Gaussian 
wave packet at various times. The top panel shows the wave packet at and just after 
the initial time. The t = lOOr plot is typical of a fractional revival where the wave 
packet reforms into a small number of recognizable features. The t = 124r is typical 
of a more random later time when the wavefunction approaches the 'flat' probability 
distribution, -Pflat('^) = V-^; indicated by the horizontal arrow. The t ~ 200r cases 
show a special fractional revival where an initially central {xq = 0) wave packet reforms 
with the initial value of Axq. Finally, the reformation of the t = wave packet at the 
half- revival time, t = T/2 = 400r is indicated, as is the reformation of the t = t/8 
packet, but flipped about the center of the well as in Eqn. Q. 

Fig. 3. Momentum-space probability distributions, t)P versus p, corresponding to the 
sample times shown in Fig. 2; the vertical dashed lines indicate the values p = ±po = 
±nQTTh/L = ±4007r. Note the reversal in momentum values at the half-revival time, 
t = T/2 = 400r. 

Fig. 4. Average value for position, {x)t (top), and spread in position, Axt (bottom), for the 
quasi-Gaussian wavepacket defined by the parameters in Eqns. (HH), (PUj) . and (^1]) 
over a time interval corresponding to the first ten classical periods. In the bottom 
figure, the first three multiples of the spreading time, to, are shown (vertical dashed 
lines). The horizontal dashed line corresponds to the initial spread, Axq = 0.05. The 
dotted curve is given by the expression for the time-dependent Axt for free-particle 
Gaussian wavepackets given by Eqn. (fTH|). 
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Fig. 5. Classical motion (x(t) and v{t) versus t) of a point particle in an infinite square well 
potential, starting at the center of the well with constant speed vq, moving in the 
positive-x direction. 

Fig. 6. Same as Fig. 4, but for the average value of momentum, {p)t (top), and spread in 
momentum, Apt (bottom). In the bottom plot, the initial momentum spread, Apo, is 
shown as a horizontal dashed line, as is the maximum value, Ap = +po, corresponding 
roughly to two equal peaks in |0(p, t)P at p = ±po as in Eqn. 

Fig. 7. Average value of position, {x)t versus t, over one revival time (T = 800r). The average 
value of position for t = {n + l/8)r (sohd) and {n + 5/8)r (dashed) are plotted in 
the top figure. (For values given hj t = nr, one obtains the horizontal solid line 
corresponding to {x)t = 0.5L.) The revival at t = T = SOOr where expectation values 
return to their t ^ values is evident, as is the half-revival at t = T/2 = 400t 
where the expectation values are mirrored about the center of the well as predicted by 
Eqn. (jH}. The value of {x)t for all times during the interval (116r, 164r) is shown on 
the bottom plot, indicating the small oscillations about the 'fiat' value of {x) = 0.5L 
during most of the collapsed phase. 

Fig. 8. Uncertainty or spread in position, Axt versus t (evaluated at times given hj t = nr) 
over one revival time (T = 800r). The 'fiat' value of Axgg^^ = L/VT2 ^ 0.288L 
is shown as the dashed horizontal line. The spread in position initially follows the 
dashed curve (indicated by Ax^^.^^) given by the free-particle expression in Eqn. (fTHjl . 
but then 'turns over' and approaches the 'fiat' value. The time interval (116r, 164r) 
is also shown to indicate the oscillations around the 'fiat' value during the collapsed 
phase. The full revival at t = T = 800r as well as the half-revival at t = T/2 = 400r 
are evident, as are the special quasi-revivals at T/4 = 200r and 3r/4 = 600r which 
appear in the xq = L/2 case only. 

Fig. 9. Average value {{p)t versus t) and uncertainty {Apt versus t) in momentum for the 
quasi-Gaussian wave packet over one revival time. The reversal in momentum seen 
in {p)t and the reformation of the wave packet to its initial width are obvious at the 
half-revival time, t = T/2 = 400r. The spread in momentum corresponding to the 
'flat' distribution of probability in position-space {Ap = +po) is also evident over a 
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large fraction of the collapsed phase. 

Fig. 10. Plot of Axt versus t showing the approach to the 'fiat' distribution for various values 
of the initial spread, Axq. In each case, the curves initially follow the free-particle 
Gaussian expression in Eqn. (fTS|) . but then 'turn over' and saturate at Ax = L/\/T2. 

Fig. 11. Plots of {x)t (top), {p)t (middle), and Apt (bottom) showing the approach to the values 
of these quantities corresponding to the 'flat' distribution. Note that the characteristic 
time for flattening is the same in all three cases and that it scales as indicated in 
Eqn. (P^. namely tgg^^ oc Axq. 
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